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Objectives: In British Columbia (BC), Canada, neuraminidase inhibitors (NIs) were publicly funded during the 2009 
A(HlNl)pdm09 pandemic for treatment of high-risk patients and/or anyone with moderate-to-severe illness. 
We assessed antiviral effectiveness (AVE) against hospitalization in that context. 

Methods: A population-based cohort study was conducted using linked administrative data. The cohort included 
all individuals living in BC during the study period (1 September to 31 December 2009) with a diagnostic code con- 
sistent with influenza or pandemic H1N1. The main study period pertained to the second-wave A(HlNl)pdm09 
circulation (1 October to 31 December 2009), with sensitivity analyses around the more specific pandemic peak 
(18 October to 7 November). Exposure was defined by same-day NI prescription. The main outcome was 
all-cause hospitalization within 14 days of the outpatient influenza diagnosis. Cox proportional hazards models 
assessed AVE with 1 : 1 propensity-score matching and covariate adjustment. 

Results: After matching, there were 304/58061 Nl-exposed and 345/58061 unexposed patients hospitalized 
during the main study period. The very young [<6 months (35.0; 95% CI 16.7-73.4)], the old [65-79 years 
(13.7; 95% CI 10.1 - 18.6)] and the very old [>80 years (38.7; 95% CI 26.6-56.5)] had the highest hospitalization 
rate per 1000 patients overall. Fully adjusted AVE against all-cause hospitalization during the main study period 
was 16% (95% CI 2%-28%), similar to the pandemic peak (15%; 95% CI -4%-30%). 

Conclusions: The use of NIs was associated with modest protection against hospitalization during the 2009 
pandemic, but appeared underutilized in affected age groups with the highest hospitalization risk. 
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Introduction 

A novel strain of influenza A(HlNl)pdm09 virus was first detected 
in Mexico 1 and the USA 2,3 in April 2009. After its initial detection, the 
pandemic virus spread to many parts of the world, including 
Canada. 4 " 6 In British Columbia (BC), Canada, the first wave of 
A(HlNl)pdm09 activity during spring/summer of 2009 was 
limited, but was followed by a second, more substantial and wide- 
spread wave in the autumn that began in early October, peaked 
during the last week of October and resolved by the end of 
2009. 7,8 Alongside other provinces of Canada, BC provided two 
pharmaceutical interventions free-of-charge as part of public 
health population efforts to mitigate overall impact of the pan- 
demic: an AS03-adjuvanted vaccine as prevention and a neur- 
aminidase inhibitor (NI) strategy as treatment. 8 

The adjuvanted A(HlNl)pdm09 vaccine used in Canada was ul- 
timately shown to be highly effective (>90%) against medically 



attended, laboratory-confirmed A(HlNl)pdm09 illness, 9 but was 
delayed in availability such that initial administration coincided with 
the pandemic peak. 8,9 Prior to the pandemic, Canada had stockpiled 
antivirals, namely the two NIs oseltamivir and zanamivir, for treat- 
ment. In anticipation of the second pandemic wave, BC released 
and distributed NIs from its emergency stockpile for the treatment 
of people at high risk of influenza illness (regardless of severity) and 
for previously healthy people experiencing moderate-to-severe 
influenza-like illness (ILI). The former could also obtain prescriptions 
in advance to be filled in the event of ILI, minimizing delay from illness 
onset to the start of treatment. Ultimately, between the weeks of 4 
October and 22 November 2009, >120000 NI prescriptions were 
filled in BC, with peak administration coinciding with other indicators 
of peak pandemic activity during the last weekof October and the first 
week of November. 8 

Although NI stockpiles are a component of pandemic prepared- 
ness in most developed countries, the 2009 pandemic was the first 
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pandemic for which widespread NI use was implemented on a 
population level with a goal of reducing serious disease and the 
associated healthcare burden. However, studies to assess antiviral 
effectiveness (AVE) against severe complications of influenza- 
seasonal or pandemic— remain largely lacking. Here, we assess 
AVE against hospitalization among patients clinically diagnosed 
with influenza during the autumn 2009 wave of A(HlNl)pdm09 
in BC, Canada. 

Methods 

Study design and population 

We conducted a retrospective cohort study using linked, administrative 
healthcare data extracted from the Medical Services Plan (MSP) billing infor- 
mation, the Hospital Separations and the PharmaNet prescription data- 
bases, provided by the BC Ministry of Health. Each eligible resident of BC is 
assigned a unique patient identifier, the personal health number, which is 
captured in all the databases and was used to link patients' records 
across the various data files. The final anonymized dataset was sent to 
the BC Centre for Disease Control in Vancouver, BC, for analysis. This study 
received approval from the University of British Columbia Research Ethics 
Board. 

The cohort(s) included all BC residents since 1 September 2009 with an 
outpatient clinical diagnosis of influenza defined by an MSP fee-service 
billing code specific for A(HlNl)pdm09 or referring to the International 
Classification of Diseases (ICD) 9th revision for influenza (ICD-9 code 
487). The date of clinical influenza diagnosis became the referent for estab- 
lishing exposure and outcome status. If the patient had more than one MSP 
diagnosis of influenza since 1 September, only the first was counted and 
used as the referent for both exposure and outcome classification. 

The main study period spanned clinical influenza diagnosis during the 
dominant second-wave A(HlNl)pdm09 activity (1 October to 31 Decem- 
ber 2009), with sensitivity analyses conducted around the more specific 
peak period (18 October to 7 November) and the broader, but less specific, 
autumn period (1 September to 31 December) commencing prior to sub- 
stantial A(HlNl)pdm09 second-wave circulation in BC. 

Antiviral exposure was defined by the filling of an NI (oseltamivirorzana- 
mivir) prescription on the same referent date (day 0), as obtained from 
PharmaNet, a population-based prescription drug database that captures 
all outpatient prescription drugs dispensed in the province, regardless of 
the payer/insurer. The referent date was chosen in the absence of informa- 
tion on actual illness onset, recognizing that delay to medical visit and 
prescription would already tend to underestimate AVE. Those filling an NI 
prescription before or after the referent date were excluded. 

The main outcome was all-cause hospitalization within 14 days of the 
referent date, obtained from the Hospital Separations database that 
includes records from all acute care inpatient visits and long-term care 
holding beds, extended care beds, rehabilitation beds and discharge plan- 
ning units managed by hospitals. In addition to all-cause hospitalization, 
we explored outcomes of pneumonia or influenza (P&I; ICD-10 codes: 
J10-18) and acute respiratory diseases (ARD; ICD-10 codes: J00-06 and 
J20-22), including where these were coded as the primary cause or any- 
where on the hospital discharge record. 

Statistical analysis 

Cox proportional hazards models assessed AVE with 1:1 propensity score 
matching and covariate adjustment. Since the assignment of subjects to 
NI treatment and non-treatment groups was not random during the pan- 
demic, confounding factors may bias treatment effects. Therefore, propen- 
sity score matching, estimated by multivariate logistic regression with the 
smallest Akaike's information criterion, was used to minimize such bias 
and balance the baseline and clinical characteristics between the two 



groups. One-to-one matching of Nl-exposed to unexposed individuals 
based on the propensity score was performed by using the 'greedy nearest- 
neighbour' algorithm. 10 Separate propensity score-matched cohorts were 
constructed for the main study period and sensitivity analyses. 

Participant profiles and hospitalization incidences were derived by NI 
exposure status before and after propensity score matching, further strati- 
fied by relevant subgroups. Multivariable Cox regression analyses compared 
the hazard of hospitalization in Nl-exposed subjects with the hazard of hos- 
pitalization in unexposed subjects in the propensity score-matched cohorts 
by the hazard ratio (HR), adjusted for relevant covariates (see below). AVE 
was derived as (1 - HR) x 100%. 

The baseline covariates used in propensity score construction and Cox 
regression analysis included age (<6 months, 6-11 months, l-4years, 
5-9 years, 10-19 years, 20-49 years, 50-64 years, 65-79 years and 
>80 years), gender, health services delivery area, number of physician 
visits on the referent date (1 and >1), number of physician visits (0, 1, 2-4 
and >5) and hospitalizations (0, 1, 2-3 and >4) within 6 months prior 
to the referent date, use of immunosuppressives within +30 days of 
the referent date, cardiorespiratory conditions (myocardial infarction, 
congestive heart failure, chronic obstructive pulmonary disease, asthma, 
cystic fibrosis and bronchopulmonary dysplasia), immunosuppressive condi- 
tions (connective tissue disease -rheumatic disease, cancer, metastatic 
carcinoma and HIV), metabolic conditions (diabetes mellitus and other 
metabolic diseases), neurological conditions (stroke and other conditions) 
and other conditions (liver disease and renal disease). Covariates were 
entered into the model via a stepwise multivariate Cox regression model 
(entry criteria of P<0.1 and staying criteria of P< 0.05) and some covariates 
were recategorized when needed. 

Information on comorbid conditions was extracted according to diag- 
nostic codes from the MSP database and Hospital Separations database 
within 2 years from the reference date (ICD-9/ICD-10 codes applied avail- 
able upon request). Adjustment for comorbidity was dichotomized as 
'yes'/'no' to any of the above chronic conditions, but also explored based 
on the Charlson index for which a score of 0 indicates no comorbidity. 11 
Patients considered immunosuppressed on the basis of therapy were 
defined by PharmaNet record of the following prescriptions: antirheumatic 
drugs, oral glucocorticoids, antirejection medication and chemotherapeu- 
tic agents. 

Statistical significance in this study was defined as P< 0.05. SAS version 
9.2 (SAS, Cary, NC, USA) was used for all statistical analyses. 

Results 

Participants 

Table SI (available as Supplementary data atJACOnline) provides a 
summary of participant profiles before propensity score matching 
and Table 1 after propensity score matching according to exposure 
and outcome status for the main analysis period of 1 October to 31 
December 2009. Table 2 compares hospitalization events by anti- 
viral exposure for the main and sensitivity analyses before and 
after propensity score matching. 

There were 227755 people who had a physician visit related to 
influenza between 1 Septemberand31 December 2009 (Figure 1). 
Of those, 10503 (4.6%) were removed due to missing information 
(n = 568) or because they filled an NI prescription after the 
referent date (n = 9935). Of the remaining 217 252 subjects, 27% 
(n = 58978, only 86 were given zanamivir) met NI exposure criteria. 
The propensity score-matched cohort for the broadest period 
(1 September to 31 December 2009) included 58775 records per 
exposure group (a total of 203 Nl-treated subjects were removed 
due to extreme propensity scores or no matching from the 
non-NI-treated pool). For the main analysis period spanning 



1398 



Table 1. Participant profile by exposure and hospitalization, main analysis period (1 October to 31 December 2009), after propensity score matching 



Hospitalized Not hospitalized Hospitalization rates per 1000 (95% CI) 

Baseline 

characteristics/ No NI NI (column %) overall No NI NI (column %) overall 

category (column %) (row%) (column %) P value (column %) (row%) (column %) P value No NI NI overall 



Age 



<6 months 


5 (1.4%) 


2 


(0.7%) (28.6%) 


7 (1.1%) 


0.088 


97 (0.2%) 


96 (0.2%) (49.7%) 


193 (0.2%) 


1.000 


49.0 (20.4, 117.8) 


20.4 (5.1,81.6) 


35.0(16.7,73.4) 


6-11 months 


1 (0.3%) 


1 


m 7o/\ /cn noc\ 
,UJ /o) \DV.V /o) 


1 in 7o/„\ 
Z (U.7 /o) 




7C7 in col\ 
70Z (u.o /oj 


7C.7 in co/\ fen 7o/„\ 

30/ (U.D 10} (7U.7 /Oj 


77Q /n co/\ 

/ zy \u.o /o) 




1 0 /n /, 1 n c\ 

z.o (u.^t, ly.oj 


7 7 in /, 1 o t.\ 


17/07 7 n n\ 

z. / (u. /, lu.yj 


1-4 years 


24 (7.0%) 


21 


Ic no/ \ It c ~ioi \ 
.6.9 /0) (46. / /o) 


/ c ic no/ \ 
45 (6.9%) 




/ one 10 / 0/ \ 
4825 (8.4 /o) 


/ 07 / to 1 0/ \ len no/ ) 
4834 (8.4%) (5U.U %) 


nccn I o / 0/ \ 
9659 (o.4 /o) 




4.9 (3.3, 7.4) 


4.3 (2.8, 6.6) 


4.6 (3.5, 6.2) 


5-9 years 


25 (7.2%) 


18 


(5.9%) (41.9%) 


43 (6.6%) 




7732 (13.4%) 


7754 (13.4%) (50.1%) 


15486 (13.4%) 




3.2 (2.2,4.8) 


2.3 (1.5, 3.7) 


2.8 (2.1, 3.7) 


10-19 years 


36 (10.4%) 


26 


(8.6%) (41.9%) 


62 (9.6%) 




12 207 (21.2%) 


12175 (21.1%) (49.9%) 


24382 (21.1%) 




2.9 (2.1,4.1) 


2.1 (1.5,3.1) 


2.5(2.0,3.3) 


20-49 years 


174 (50.4%) 


138 


(45.4%) (44.2%) 


312 (48.1%) 




24531 (42.5%) 


24 587 (42.6%) (50.1%) 


49118 (42.5%) 




7.0 (6.1,8.2) 


5.6 (4.7, 6.6) 


6.3(5.6, 7.1) 


50-64 years 


57 (16.5%) 


52 


(17.1%) (47.7%) 


109 (16.8%) 




6111 (10.6%) 


6108 (10.6%) (50.0%) 


12 219 (10.6%) 




9.2 (7.1, 12.0) 


8.4(6.4,11.1) 


8.8(7.3,10.7) 


65-79 years 


14 (4.1%) 


28 


(9.2%) (66.7%) 


42 (6.5%) 




1510 (2.6%) 


1507 (2.6%) (50.0%) 


3017 (2.6%) 




9.2(5.4,15.5) 


18.2 (12.6, 26.4) 


13.7(10.1,18.6) 


>80 years 


9 (2.6%) 


18 


(5.9%) (66.7%) 


27 (4.2%) 




341 (0.6%) 


329 (0.6%) (49.1%) 


670 (0.6%) 




25.7 (13.4, 49.4) 


51.9 (32.7, 82.3) 


38.7 (26.6,56.5) 


Sex 
female 


215 (62.3%) 


181 


(59.5%) (45.7%) 


396 (61.0%) 


0.469 


31541 (54.6%) 


31612 (54.7%) (50.1%) 


63153 (54.7%) 


0.774 


6.8 (5.9, 7.7) 


5.7 (4.9, 6.6) 


6.2(5.6, 6.9) 


male 


130 (37.7%) 


123 


(40.5%) (48.6%) 


253 (39.0%) 




26175 (45.4%) 


26145(45.3%) (50.0%) 


52320(45.3%) 




4.9 (4.2, 5.9) 


4.7 (3.9, 5.6) 


4.8 (4.3, 5.4) 


Immunosuppressive drug use 
























yes 


39 (11.3%) 


50 


(16.4%) (56.2%) 


89 (13.7%) 


0.057 


1485 (2.6%) 


1520 (2.6%) (50.6%) 


3005 (2.6%) 


0.531 


25.6 (18.7,35.0) 


31.8 (24.1,42.0) 


28.8 (23.4,35.4) 


no 


306 (88.7%) 


254 


(83.6%) (45.4%) 


560 (86.3%) 




56231 (97.4%) 


56237(97.4%) (50.0%) 


112468 (97.4%) 




5.4 (4.8, 6.1) 


4.5 (4.0, 5.1) 


5.0 (4.6, 5.4) 


No. of same-day GP visits 
























1 


309 (89.6%) 


289 


/or 10/ \ 1 1 o 10/ \ 
,95.170) [HO. 370) 


crno /m 1 0/ \ 

byy (y2.i%) 


0.009 


c~? / a /no co/ \ 

5/452 (yy.5%) 


chi cc /nn co/ \ icr\ no/ \ 
5/455 (yy.5%) pU.Uvo) 


11/ nm /nn co/ \ 

iiM-yu/ (yy.5%) 


0.111 


c 0 1/ 0 c n\ 
5.3 [h.o, b.\J) 


5.0 (4.5, 5.6) 


5.2 (4.8, 5.6) 


2 


32 (9.3%) 


12 


(3.9%) (27.3%) 


44 (6.8%) 




263 (0.5%) 


298 (0.5%) (53.1%) 


561 (0.5%) 




108.5 (76.7, 153.4) 


38.7 (22.0, 68.2) 


72.7 (54.1,97.7) 


>3 


4 (1.2%) 


3 


(1.0%) (42.9%) 


7 (1.1%) 




1 (0.0%) 


4 (0.0%) (80.0%) 


5 (0.0%) 




800.0 (300.2,2131.6) 


428.6 (138.2,1328.8) 


583.3 (278.1, 1223.6) 


No. of past GP visits 


























0 


321 (93.0%) 


2 9 b 


/Q7 no/ \ 1 /. 1 no/ \ 
,y /.U/o; (4/.y7oJ 


tit fa/, no/ \ 


0.095 


CCQ~!C IQQ 10/ \ 

bo y / d (yo./ /oj 


C~! C\C\ fQQQO/WCn no/ >, 
D/Lol l,yo.o /oj PU.U /oj 


11/. mc IQQ 70/ \ 


0.528 


c c ic n c ~2\ 


r 1 //. c c. Q\ 

D.L [H.O, D.O) 


r /, ic C\ C Q\ 
D.4 p.U, D.O) 


1 


17 (4.9%) 


5 


/1 co/ w*n "70/ \ 
,1.0%) h) 


22 (3.4%) 




cm /1 no/ \ 
byz (l.U7o) 


c cr\ /1 no/ \ // n no/ \ 
56y (l.UTo) (4y.(J7o) 


1 1 c 1 /1 nn/ \ 
1 161 (l.U7o) 




-)"7n/i"7 / / / n\ 
2 /.y (1 /A, 44.9) 


0 1 /1 c ">n n\ 
0./ (3.b, 2U.y) 


lo.b (12.2, 2o.2) 


2-4 


5 (1.4%) 


2 


/A "70/ \ /TO CO/ \ 

.U. /7o) (2o.b7o) 


"7/1 1 0/ \ 
1 (1.1%) 




1 ~i 1 /n in/ 1 
12 / (U.27o) 


1 1 n M in/ \ / / c / o/ \ 
llU (U.27o) (4oA7o] 


~t "7 /n m/ \ 
23 / (U.27o) 




■3"7n/ico m n\ 
3/.y (15. 0, 91.UJ 


1 1 n // c "71 / \ 
1 /.9 (4.5, / 1.4) 


1 O *7 M 13 *7 Cr\ 1\ 

lo.i (13./, oU.2) 


>5 


2 (0.6%) 


2 


(0.7%) (50.0%) 


4 (0.6%) 




22 (0.0%) 


27 (0.0%) (55.1%) 


49 (0.0%) 




83.3 (20.8, 333.2) 


69.0 (17.2, 275.8) 


75.5 (28.3, 201.1) 


No. of past hospitalizations 
























0 


271 (78.6%) 


2 59 


/70 co/ \ I / c no/ i 
JQXrfa) |46. 9 /o) 


c i n t~/o co/ i 
51U (/o.b/oj 


0.632 


c / c 07 mi enn 
54 581 (94.6/o) 


c / c t n In/ co/ \ fen on/ i 
54 64U (94.6%) (50.0 /o) 


-inn 77-1 /n/ co/ \ 
10y 221 (94.6%) 


0.758 


4.9 (4.4, 5.6) 


/ / 11 0 / n\ 
4.4 (3.8, 4.9) 


4.6 (4.3, 5.1) 


1 


47 (13.6%) 


39 


/ T 7 oo/ \ lie 70/^ 
,12.070) (45.3%) 


oc 11 1 00/ 1 
86 (13.3%) 




7 77 1 // "70/ I 

2/11 (4./%) 


7cm // 70/ i // n 00/ \ 
2692 (4./ 70) (4y.8 /o) 


c / no // 70/ \ 
5403 (4. / /o) 




1 /.U (12. 0, 22. /) 


i / 7 '1 n / in 
14.3 (1U.4, 19.5) 


15./ (12./, 19.4) 


2-3 


19 (5.5%) 


22 


/7 10/ \ 10 -70/ \ 
J. 1 /o) (53./ /o) 


/ i /C 70/ \ 

41 (6.3/o) 




7 m /n "70/ i 
39/ (0. / 70) 


7nn /n 70/ w/ n co/ 1 
39U (0. /70) (49.6 /oJ 


707 /n 70/ \ 

/8/ (U./70) 




/ c i nn 1 "71 c\ 
45./ (29.1, / l.o) 


C7/ /7C7 Ol 1\ 

53.4 (35.2, 81.1) 


/ n c i~} c c c~i i\ 
4y.5 (36.5, fa/. 2) 


>4 


8 (2.3%) 


4 


/i in/ 1 i~> ~> 00/ \ 
.1.3 /o) (33.3 /o) 


77/1 00/ \ 
12 (I.5 /0) 




m /n no/ 1 
2 / (0.0 70) 


7 c /n 7 0/ \ /cc co/ 1 
35 (O.l/o) (56.5 70) 


a In 1 0/ \ 
62 (U.l /o) 




22o.6 (114.3, 45/. 1) 


im c no c 777 7^ 
102.6 (38.5, 2/3.3) 


1 c n t /m 1 ioc c\ 
162.2 (y2.1, 2o5.5) 


Charlson index 


























0 


247 (71.6%) 


230 


(75.7%) (48.2%) 


477 (73.5%) 


0.311 


54275 (94.0%) 


54244 (93.9%) (50.0%) 


108519(94.0%) 


0.195 


4.5 (4.0, 5.1) 


4.2 (3.7, 4.8) 


4.4 (4.0, 4.8) 


1-2 


76 (22.0%) 


59 


(19.4%) (43.7%) 


135 (20.8%) 




3155 (5.5%) 


3195 (5.5%) (50.3%) 


6350 (5.5%) 




23.5(18.8, 29.5) 


18.1 (14.0,23.4) 


20.8(17.6,24.6) 


3-5 


13 (3.8%) 


13 


(4.3%) (50.0%) 


26 (4.0%) 




213 (0.4%) 


252 (0.4%) (54.2%) 


465 (0.4%) 




57.5(33.4, 99.1) 


49.1 (28.5, 84.5) 


53.0 (36.1,77.8) 


6-9 


8 (2.3%) 


2 


(0.7%) (20.0%) 


10 (1.5%) 




72 (0.1%) 


62 (0.1%) (46.3%) 


134 (0.1%) 




100.0 (50.0, 200.0) 


31.3 (7.8, 125.0) 


69.4 (37.4,129.1) 


>10 


1 (0.3%) 


0 


(0.0%) (0.0%) 


1 (0.2%) 




1 (0.0%) 


4 (0.0%) (80.0%) 


5 (0.0%) 




500.0 (70.4,3549.7) 


0 


166.7 (23.5,1183.2) 


Age/Charlson index 


























<6 months/>l 


5 (1.4%) 


2 


(0.7%) (28.6%) 


7 (1.1%) 


<0.001 


97 (0.2%) 


95 (0.2%) (49.5%) 


192 (0.2%) 


<0.001 


49.0 (20.4,117.8) 


20.6 (5.2,82.4) 


35.2(16.8, 73.8) 


6-11 months/0 


0 (0.0%) 


0 


(0.0%) (0.0%) 


0 (0.0%) 




0 (0.0%) 


1 (0.0%) (100.0%) 


1 (0.0%) 










<6 months/0 


1 (0.3%) 


1 


(0.3%) (50.0%) 


2 (0.3%) 




354 (0.6%) 


362 (0.6%) (50.6%) 


716 (0.6%) 




2.8 (0.4,20.0) 


2.8 (0.4, 19.6) 


2.8 (0.7,11.1) 


6-11 months/>l 


0 (0.0%) 


0 


(0.0%) (0.0%) 


0 (0.0%) 




8 (0.0%) 


5 (0.0%) (38.5%) 


13 (0.0%) 










1-4 years/0 


19 (5.5%) 


17 


(5.6%) (47.2%) 


36 (5.5%) 




4515 (7.8%) 


4515(7.8%) (50.0%) 


9030 (7.8%) 




4.2 (2.7,6.6) 


3.8 (2.3,6.0) 


4.0 (2.9,5.5) 


l-4years/>l 


5 (1.4%) 


4 


(1.3%) (44.4%) 


9 (1.4%) 




310(0.5%) 


319 (0.6%) (50.7%) 


629 (0.5%) 




15.9 (6.6,38.1) 


12.4 (4.6,33.0) 


14.1 (7.3,27.1) 



Continued 
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1 October to 31 December 2009, there were 58061 per group; for 
the most specific pandemic peak period spanning 18 October to 
7 November, there were 36771 records per group. 

Before propensity score matching, the two cohorts (NI exposed 
and unexposed) differed significantly on almost all baseline char- 
acteristics (11 showed a Pvalue of <0. 0002 and only cardiorespira- 
tory condition had a P value of 0.7). These variables were used to 
derive propensity scores upon which the treatment groups were in- 
dividually matched in each analysis period. After 1:1 propensity 
score matching, no baseline characteristics, including those con- 
sidered as possible confounders, showed significant differences 
between groups. The distribution of all baseline covariates 
was completely balanced between Nl-exposed and unexposed 
groups by the propensity score matching. Since only 203 subjects 
(0.03%) from the Nl-treated subjects were lost during the match- 
ing algorithm, the final matching sample retains the representa- 
tiveness of the population. Both before and after propensity score 
matching, data showed similar patterns in Nl-exposed and 
unexposed groups with respect to the distribution of intervals 
between influenza visit and subsequent hospitalization. More 
than 50% of hospitalized subjects were admitted by day 3. 

Overall and among subjects in both Nl-exposed and unexposed 
groups, the highest hospitalization rates after propensity score 
matching were in the very young (<6 months old) as well as the 
old (65-79 years old) and the very old (>80 years old) (Table 1). 
Overall rates of hospitalization per 1000 patients in the exposed 
and unexposed cohorts, within 2 weeks of an outpatient influenza 
diagnosis, were significantly higher in these age groups than in any 
other: 35.0 (95% CI 16.7-73.4), 13.7 (95% CI 10.1-18.6) and 38.7 
(95% CI 26.6-56.5), respectively (Table 1). These ages comprised 
0.2%, 2.6% and 0.6% of participants with outpatient influenza 
diagnosis. 

About 6% of subjects with an outpatient physician diagnosis of 
influenza who were not subsequently hospitalized (i.e. within 
14 days) had an underlying comorbidity (Table 1). Conversely, 
among hospitalized patients, about one-quarter of the exposed 
and unexposed groups had an underlying comorbidity and this 
was mainly due to age groups 20-49 and 50-64 years. 
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AVE 

The crude and adjusted baseline HRs and the 95% CIs for all-cause 
hospitalization associated with the use of antivirals are illustrated 
in Table 3. During the main analysis period (Table 4), spanning 1 
October to 31 December 2009, antivirals were associated with a 
statistically significant reduction of 16.1% in the riskof hospitaliza- 
tion (HR 0.839; 95% CI 0.719-0.980), comparable to the narrower, 
but more specific, pandemic peak period of 18 October to 7 Novem- 
ber 2009 (AVE 15%; 95% CI -4%-30%). For the broader, but less 
specific, period spanning 1 September to 31 December 2009, AVE 
was higher (36%; 95% CI 20%-49%), but there was greater vari- 
ability in this estimate. During the main analysis period, the use 
of more specific causes of hospitalization (P&I or ARD) also para- 
doxically resulted in lower AVE estimates. 

For the main analysis period, crude AVE among those with 
comorbidity was significant at 52% (95% CI 29% -68%) as was 
the fully adjusted model (59%; 95% CI 39%-73%). For those 
individuals without comorbidity, the crude and the fully adjusted 
AVE were non-significant at 6% (95% CI -14% -20%) and 6% 
(95% CI -13%-21%), respectively. 



Table 2. Crude all-cause hospitalization and death rate per 100000 person-days by antiviral exposure status and analysis period before and after propensity score matching 











Before propensity score matching 




After propensity score matching 






Antiviral 




no. of 








no. of 






Time period 


Outcome 


exposure 


n 


events 


person-days 


rate (95% CI) 


n 


events 


person-days 


rate (95% CI) 


Main period (1 October 


hospitalization 


total 


213022 


1343 


2 969 756 


45.22 (42.87, 47.71) 


116122 


649 


1619 569 


40.07(37.10, 43.28) 


to 31 December) 




antiviral 


58271 


310 


812 945 


38.13 (34.12,42.62) 


58 061 


304 


810 087 


37.53 (33.54,41.99) 






no antiviral 


154 751 


1033 


2 156811 


47.89 (45.06, 50.91) 


58 061 


345 


809 482 


42.62 (38.35,47.36) 




death 


total 


213 022 


124 


6 388 512 


1.94 (1.63,2.31) 


116 122 


28 


3 483 160 


0.80 (0.56,1.16) 






antiviral 


58271 


11 


1 747 938 


0.63 (0.35, 1.14) 


58 061 


11 


1 741 638 


0.63 (0.35, 1.14) 






no antiviral 


154 751 


113 


4 640 574 


2.44 (2.03,2.93) 


58 061 


17 


1 741 522 


0.98 (0.61,1.57) 


Sensitivity analysis 






















peak period (18 October 


hospitalization 


total 


115037 


724 


1 603 762 


45.14(41.97, 48.55) 


73 542 


380 


1 025 985 


37.04 (33.49, 40.96) 


to 7 November) 




antiviral 


36 962 


184 


515 755 


35.68 (30.88, 41.22) 


36 771 


175 


513188 


34.10(29.40,39.55) 






no antiviral 


78 075 


540 


1088 007 


49.63 (45.62, 54.00) 


36 771 


205 


512 797 


39.98 (34.86,45.84) 




death 


total 


115037 


60 


3450 017 


1.74 (1.35,2.24) 


73 542 


15 


2 205 972 


0.68 (0.41,1.13) 






antiviral 


36 962 


7 


1 108 726 


0.63 (0.30, 1.32) 


36 771 


7 


1 102 996 


0.63 (0.30, 1.33) 






no antiviral 


78 075 


53 


2 341 291 


2.26 (1.73, 2.96) 


36 771 


8 


1 102 976 


0.73 (0.36, 1.45) 


full autumn (1 September 


hospitalization 


total 


217252 


1398 


3 028415 


46.16(43.81,48.65) 


117 550 


668 


1639 264 


40.75 (37.77, 43.96) 


to 31 December) 




antiviral 


58 978 


317 


822 790 


38.53 (34.51,43.01) 


58 775 


310 


820 039 


37.80 (33.82,42.25) 






no antiviral 


158 274 


1081 


2 205 625 


49.01 (46.17, 52.02) 


58 775 


358 


819225 


43.70 (39.40, 48.47) 




death 


total 


217252 


134 


6 515 214 


2.06 (1.74, 2.44) 


117550 


31 


3 525 974 


0.88 (0.62,1.25) 






antiviral 


58 978 


11 


1 769148 


0.62 (0.34, 1.12) 


58 775 


11 


1 763 058 


0.62 (0.35, 1.13) 






no antiviral 


158 274 


123 


4 746 066 


2.59 (2.17, 3.09) 


58 775 


20 


1 762 916 


1.13 (0.73, 1.76) 
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Cohort n = 227 755 with MSP visit related to influenza between 
1 September and 31 December 2009 



Subjects were removed: 568 subjects due to 
unknown demographics and 9935 subjects due 
to antiviral received after the index date 



n=217252 



NI treated 
n=58978 



NI untreated 
n=158274 



Propensity scores were created using logistic regression 
Non-overlapping scores between Nl-treated and untreated 
groups were removed before 1:1 matching 



NI treated 
n=58974 



NI untreated 
n=158265 



Perform 1:1 propensity scores matching 
(199 Nl-treated subjects were lost due to no matching 
from untreated groups) 



NI treated 
n=58775 



NI untreated 
n=58775 



117 550 subjects were used in 
proportional HR analysis 



Figure 1. Propensity score matching cohort flow chart during second wave. 



Discussion 

Our large, population-based study found that antiviral use, as mea- 
sured by receipt of a prescription for oseltamivir or zanamivir after 
a physician diagnosis of influenza or A(HlNl)pdm09 pandemic, 
was associated with a 15%- 36% decrease in all-cause hospital- 
ization. Our study is unique in that we used a population-based 
cohort to evaluate the effectiveness of antivirals; few studies 
have evaluated AVE in an entire cohort 12,13 and most used case 
series. 14 " 18 Nevertheless, our results are consistent with other 
findings that show that NIs are useful during pandemics when 



taken within the recommended 48 h. These results support the 
use of antiviral medications as a key mitigation strategy. 

Although the A(HlNl)pdm09 pandemic virus was more likely to 
infect the younger population, previous studies have showed that 
when patients aged >65 years were infected with influenza, they 
were more likely to be admitted into hospital and die from it. 19 Our 
findings were similar and we saw the highest hospitalization rates 
in the very young (i.e. <6 months of age) and in older adults 
(i.e. >65 years of age). About 6% of subjects with an outpatient 
physician diagnosis of influenza who were not subsequently hospita- 
lized (i.e. within 14 days) had an underlying comorbidity (Table 1). 
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Table 3. Crude and adjusted HR for all-cause hospitalization, main 
analysis period (1 October to 31 December 2009) 





HR (95% CI) 


P value 


Crude estimate 


0.881 (0.755, 1.027) 


0.1061 


Individual covariate adjustment 






nnp nmi in 

y ui uuu 


0 881 (0 755 1 0281 


0.1067 


gender 


0.880 (0.755,1.027) 


0.1053 


health authority 


0.880 (0.754,1.027) 


0.1047 


immunosuppressive agent use 


0.877 (0.752, 1.024) 


0.0964 


number of past GP visits 


0.883 (0.756, 1.030) 


0.1123 


number of past hospitalizations 


0.881 (0.755, 1.028) 


0.1071 


number of same-day GP visits 


0.738(0.590,0.923)° 


0.0078 


Charlson index 


0.743 (0.594,0.929)° 


0.0093 


vaccine availability 13 


0.864 (0.740, 1.009) 


0.0643 


Fully adjusted estimate 


0.839 (0.719,0.980) 


0.0270 



Sample size (no. of hospitalizations): antiviral, 58061 (304); non-antiviral, 
58061 (345). 

indicating model was adjusted to non-proportionality, estimates unstable, 
indicating whether the influenza referent date was during the period that 
vaccine became available on 26 October 2009. 



Conversely, among hospitalized patients, about one-quarter of the 
treated and untreated groups had an underlying comorbidity. This 
high proportion was mainly contributed by adults, particularly 
those in age groups 20-49 and 50-64 years. 

As mentioned previously, a number of studies, mostly case 
series, have evaluated the benefit of antivirals in reducing the 
severity of infections (reduction in critical care admission), hospita- 
lizations and mortality during the 2009 pandemic. 12,14 " 18,20,21 
Of note, the majority of these studies have been conducted in preg- 
nant women. 15 " 18,20,21 Like our study, they all demonstrated that 
NIs decreased the risk of severe disease and hospitalization. Our 
study databases could not identify pregnant women as a 
separate category as there was a lack of appropriate ICD-9 or -10 
diagnostic codes to identify pregnancy during the pandemic. It 
would have been interesting to see the effect of the antivirals on 
this population, given that they are at higher risk of hospitalization 
and severe disease from influenza. 15,20,21 A few of our patients did 
go to ICU, but the sample size was not large enough to do propen- 
sity score matching for the exposed and unexposed groups; thus, 
the data are not presented. 

Early on in the 2009 pandemic, BC developed a targeted pro- 
gramme to encourage antiviral use in people with severe ILI 
disease and/or underlying chronic conditions. In our study, we 
saw that these patients were >30% more likely to receive antivir- 
als. The provincial stockpile of NIs was also pre-distributed across 
the province to pharmacies and remote communities, so that 
people at risk could have access to the medications in a timely 
way. This study supports the benefits of this targeted approach 
to those at highest risk and the strategy of antiviral pre-distribution. 
In mitigating the pandemic impact, the avoidance of between 16% 
and 36% of hospitalizations attributable to the antiviral strategy 
was likely a key factor in the ability of the healthcare system to 
cope with the peak of infections between October and November 
2009, before vaccine was available. 



There are several limitations to this study. First, the measure of 
antiviral exposure we used was a prescription for antivirals being 
dispensed to an individual on day 0 with the assumption that the 
patients had symptoms for <48 h prior to their physician visit. 
During our vaccine effectiveness analysis, 9 the median time from 
ILI onset to physician visit was 2-3 days and therefore we may 
have included patients who had had symptoms for >48 h prior 
to their physician visit; however, if this was the case, AVE would 
be higher than we observed. Further, we cannot be certain that 
the medication was actually taken by the recipient or that it was 
taken for a complete course. We are assuming that all the 
A(HlNl)pdm09 pandemic viruses were NI susceptible, but it 
should be noted that oseltamivir resistance is present in many 
countries and, as of 26 January 2011, 340 instances of oseltamivir 
resistance have been reported by the WHO Global Influenza Sur- 
veillance Network. 22,23 

Individuals may have received antivirals from sources other 
than a pharmacy and these would not be recorded in the Pharma- 
Net database. This was a potential primarily in remote communi- 
ties where antivirals were pre-positioned within the community 
to be dispensed by the local healthcare worker. In addition, 
PharmaNet would not have recorded antiviral use for inmates in 
federal penitentiaries. Additionally, some people would have 
received a prescription from their physician in advance of the pan- 
demic to be taken should they develop ILI. This was part of a cam- 
paign to ensure people at risk had a plan for assessment (often by 
phone) and treatment should they develop influenza during the 
pandemic. In either case, failure to capture these individuals 
would lead to conservative bias and an underassessment of AVE. 
Oseltamivir may have been administered to the 'healthier' popula- 
tion; a proxy for this could have been past vaccinations. Unfortu- 
nately, our merged datasets did not include past vaccinations as 
we do not have a complete immunization registry for adults. 
However, we did adjust for the 'healthy' adults by the number of 
past GP visits as well as the number of past hospitalizations. 

Miscoding is always a possibility when using administrative data, 
especially when related to coding for physician office visits or hospi- 
talization. This could have resulted in an over- or underestimation of 
AVE. This is particularly true for the varying estimates of AVE seen 
when we looked at more influenza-specific causes of hospitaliza- 
tions, such as acute respiratory disease and/or pneumonia. Hospital 
emergency departments in BC use both physicians who bill for their 
services and those who are paid a flat fee. As such, if only an out- 
patient visit was required and the patient was seen by a non-billing 
physician, we would have not included them in our dataset. 

Prior to the pandemic, randomized controlled trials evaluating 
the efficacy of antivirals for treatment of seasonal influenza A or 
B had shown them to shorten the course of illness when adminis- 
tered within 48 h of onset of illness, 24 " 30 although the benefit was 
greatest when treatment was initiated with 12 h of the onset of 
symptoms. 31 The use of antivirals in high-risk patients, defined as 
those over the age of 65 years or with chronic medical conditions, 
showed that they reduced the time to alleviation of symptoms by 
~0.5 - 1 day. 32 Although these studies showed a modest decrease 
in symptom duration, not many studies had evaluated antiviral 
efficacy in preventing hospitalization or mortality. A pooled analysis 
of 10 randomized controlled trials of oseltamivir used in adults with 
acute influenza showed that its use was associated with a 50% 
decline in the hospitalization rate or lower respiratory infections 
and antibiotic use declined by 26%. 33 
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Table 4. Main and sensitivity analyses of antiviral effectiveness by outcome 



All-cause hospitalizations 
n (no. of hospitalizations) AV 
n (no. of hospitalizations) non-AV 
HR estimates crude 
fully adjusted 

Primary due to P&I 
n (no. of hospitalizations) AV 
n (no. of hospitalizations) non-AV 
HR estimates crude 
fully adjusted 

Due to P & I (anywhere on discharge sheet) 
n (no. of hospitalizations) AV 
n (no. of hospitalizations) non-AV 
HR estimates crude 
fully adjusted 

Primary due to acute respiratory diseases 
n (no. of hospitalizations) AV 
n (no. of hospitalizations) non-AV 
HR estimates crude 
fully adjusted 



1 October to 31 December 2009 
(main period) 

58061 (304) 
58061 (345) 

0.881 (0.755, 1.027) (0.11) 
0.839 (0.719,0.980) (0.03) 

58061 (57) 
58061 (40) 

1.425 (0.951, 2.135) (0.09) 
1.421(0.947, 2. 131) (0.09) 

58061 (77) 
58061 (71) 

1.084 (0.785,1.497) (0.62) 
1.041 (0.753, 1.440) (0.81) 

58061 (11) 
58061 (11) 

1.000 (0.434, 2.307) (1.00) 
0.975(0.422,2.249) (0.95) 



Due to acute respiratory diseases (anywhere on discharge sheet) 
n (no. of hospitalizations) AV 58061 (18) 

n (no. of hospitalizations) non-AV 58061 (18) 

HR estimates crude 1.385 (0.678, 2.826) (0.37) 

fully adjusted 1.350 (0.661, 2.756) (0.41) 

Primary due to P & I and acute respiratory diseases 
n (no. of hospitalizations) AV 58061 (67) 

n (no. of hospitalizations) non-AV 58061 (51) 

HR estimates crude 1.314 (0.913, 1.891) (0.14) 

fully adjusted 1.257 (0.873, 1.812) (0.22) 

Due to P & I and acute respiratory diseases (anywhere on discharge sheet) 
n (no. of hospitalizations) AV 58061 (91) 

n (no. of hospitalizations) non-AV 58061 (84) 

HR estimates crude 0.816 (0.560, 1.189)° (0.29) 

fully adjusted 0.772 (0.529, 1.127)° (0.18) 



HR(95%CI) (P value) 

1 September to 31 December 2009 
(second wave) 



58775 (310) 
58775 (358) 

0.674 (0.542, 0.839) a (<0.01) 
0.639 (0.513, 0.796) a (<0.01) 

58775 (59) 
58775 (47) 

1.255 (0.856, 1.841) (0.24) 
1.221 (0.832, 1.793) (0.31) 

58775 (79) 
58775 (81) 

0.975 (0.715, 1.329) (0.87) 
0.935 (0.685, 1.275) (0.67) 

58775 (12) 
58775 (11) 

1.091 (0.481, 2.472) (0.84) 
1.078 (0.476, 2.444) (0.86) 

58775 (19) 
58775 (15) 

1.267 (0.644, 2.493) (0.49) 
1.239(0.629,2.440) (0.54) 

58775 (70) 
58775 (58) 

0.887 (0.564,1.396)° (0.61) 
0.845 (0.536,1.331)° (0.47) 

58775 (94) 
58775 (96) 

0.734(0.505,1.069)° (0.11) 
0.690 (0.474, 1.006)° (0.054) 



18 October to 7 November 2009 
(peak period) 



36771 (175) 
36771 (205) 

0.853 (0.697, 1.044) (0.12) 
0.850 (0.695, 1.041) (0.12) 

36771 (34) 
36771 (32) 

1.062 (0.656, 1.722) (0.81) 
1.035 (0.638, 1.679) (0.89) 

36771 (47) 
36771 (56) 

0.839 (0.570, 1.237) (0.38) 
0.808 (0.548,1.191)(0.28) 

36771 (6) 
36771 (4) 

1.500 (0.423, 5.315) (0.53) 
1.475 (0.416, 5.228) (0.55) 

36771 (9) 
36771 (8) 

1.125 (0.434, 2.915) (0.81) 
1.130 (0.436, 2.929) (0.80) 

36771 (39) 
36771 (36) 

1.083 (0.689, 1.704) (0.73) 
1.076 (0.684, 1.693) (0.75) 

36771 (54) 
36771 (64) 

0.844 (0.587, 1.212) (0.36) 
0.807 (0.562,1.160) (0.25) 



AV, antiviral. 

indicating model was adjusted to non-proportionality, estimates unstable. 



Numerous studies have now been published on the 
A(HlNl)pdm09 pandemic, but most of these looked at risk factors 
associated with pHlNl. These studies showed that patients at 
increased risk of hospitalization 6,19,34 " 37 and severe disease from 
pHlNl were those with underlying medical conditions (especially 
asthma and chronic obstructive pulmonary disease), 5,14,19,32 " 35 
children <2 years of age (especially those with asthma and 
neurological conditions), 14,19,38,39 obese patients (BMI >35) , 40 " 42 
pregnant women 15,20,21 and aboriginal peoples. 6,39,43 



The adjuvanted A(HlNl)pdm09 vaccine used in Canada was 
ultimately shown to be highly effective (>90%) against medically 
attended, laboratory-confirmed A(HlNl)pdm09 illness, 9 but was 
delayed in availability such that initial administration coincided 
with the pandemic peak in BC. Further, the vaccine was initially 
available in only limited amounts, requiring sequenced rollout 
beginning in the last week of October for persons with comorbidity 
<65 years of age, pregnant women and remote community resi- 
dents, followed by children <5 years of age, healthcare workers 
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and other caregivers <65 years in eariy November, then older chil- 
dren and first responders and, finally, all other BC residents begin- 
ning mid-late November. Uptake of the vaccine was ~35%-45% 
in the province overall and was highest in the elderly (age 
>65 years) in whom vaccine administration had been most 
delayed. 8 As such, although the use of the NIs may not have 
slowed down transmission of A(HlNl)pdm09, they were effective 
in reducing hospitalizations and therefore decreasing the burden 
on the acute care sector. 

Conclusions 

Antiviral use in people with influenza who were at risk of severe 
disease or complications from influenza infection was associated 
with a reduced risk of hospitalization during the 2009 influenza 
pandemic. Antiviral strategies should continue to be incorporated 
into pandemic planning for future influenza pandemics. 
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